Local Magnetic Turbulence and TeV— PeV Cosmic Ray Anisotropics 
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In the energy range from ~ 10^^ eV to ~ 10^^ eV, the Galactic cosmic ray flux has anisotropics 
both on large scales, with an amplitude of the order of 0.1%, and on scales between ~ 10° and 
~ 30°, with amplitudes smaller by a factor of a few. With a diffusion coefficient inferred from 
Galactic cosmic ray chemical abundances, the diffusion approximation predicts a dipolar anisotropy 
of comparable size, but does not explain the smaller scale anisotropics. We demonstrate here that 
energy dependent smaller scale anisotropics naturally arise from the local concrete realization of the 
turbulent magnetic field within the cosmic ray scattering length. We show how such anisotropics 
could be calculated if the magnetic field structure within a few tens of parsecs from Earth were 
known. 



Introduction. — In the TeV-PeV energy range, statisti- 
cally significant anisotropics in the distribution of Galac- 
tic Cosmic Ray (CR) arrival directions on the sky have 
been reported both on large and small scales by sev- 
eral experiments, such as Super-Kamiokande Tibet- 
Ill 01, Milagro 0,13, ARGO-YBJ and IceCube 
On the largest scales these anisotropics have an ampli- 
tude of the order of ~ 0.1%, and are smaller by a fac- 
tor of a few to ten on scales between ~ 10° and ~ 30°. 
They are detected up to £; ~ 400 TeV by IceCube 
and EAS-TOP [|. While they appear strongly energy- 
dependent 0, with Milagro even finding a localized ex- 
cess with a different spectrum Q than for the rest of the 
sky, they seem to be stable in time [lo| . 

The large scale anisotropy can be explained within the 
diffusion approximation. The transport of charged CRs 
in the Galactic magnetic field is diffusive at least up to 
E ~ 10"'^^""'^'' eV. The gyroradius of particles of momen- 
tum p and charge Ze in a magnetic field of strength B, 
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is indeed smaller than the largest length scales, Lmax = 
100 — 300 pc, on which the turbulent component of the 
Galactic magnetic field is believed to vary. CRs pre- 
dominantly scatter on those turbulent field modes whose 
wave vectors k satisfy |k| ~ 27r/rg(p). As a result, the 
CR propagation in the Galaxy resembles a random walk. 
It can be considered as a Markovian process in which 
CRs "lose" memory of their past trajectory on distance 
scales of the order of the scattering length on magnetic 
field inhomogeneities X{p), which here is essentially the 
length scale beyond which the CR propagation direction 
becomes uncorrelated with its original direction. 

In the diffusion approximation, an inhomogeneous den- 
sity n{r,p) of CRs with momenta p leads to a dipole vec- 
tor given by [ni[ii 

3 j _ 3D{p) Vn 
cq n Co n ' 

where j(r,p) = —D{p)'Vn is the CR current correspond- 
ing to the diffusion coefficient D{p) = A(p)/3, here as- 
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sumed to be homogeneous, and cq is the speed of light. 
A nearby source would lead to a dipole oriented towards 
the source. The contribution of several recent nearby 
sources to the CR fiux at Earth would be a superposi- 
tion of dipoles and thus again a dipole. Eq. ^ yields the 
estimate \S{p)\ ^ D{p)/R where R is of the order of the 
distance to the closest sources. D{p) can be inferred from 
the chemical abundance ratios of secondary to primary 
nuclei, such as boron to carbon, and from the anti-proton 
fraction This gives fits of the form [li,[ll| 
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where zo ^ 1 kpc is the scale height of the Galactic disk, 
and 5 ~ 0.45. Theoretically, (5 = 2 — a, for a turbulent 
field with spectral index a: (5=1/3 for a Kolmogorov 
spectrum. With Eq. ([2]) and zq this yields 



\S{p)\ ~ few X 10" 
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More detailed treatments such as those in Refs. |16l. Il7| 
give similar values and have found that such large scale 
anisotropies contain the signature of the few most nearby 
and recent sources. For a given source distribution, 
may strongly fiuctuate around its average scaling with 
[l7| . This can explain why the anisotropies mea- 
sured by IceCube 0] are smaller at 400 TeV than at 
20 TeV. Other works have proposed that the large scale 
anisotropies may be due to a combined effect of the reg- 
ular and turbulent Galactic Magnetic Field (GMF) jlSj . 
or to local uni- and bi-directional infiows [l9| . 

However, the intermediate and small scale anisotropies 
remain hard to explain. The nearest objects believed to 
contribute to Galactic CRs (such as Vela, Gcminga or the 
Gum nebula) are at distances > 200 pc, and thus many 
gyroradii away, sec Eq. ([1]). Therefore, anisotropics at 
scales smaller than the dipole should not correlate with 
the directions towards possible sources. This would only 
be the case in the presence of additional effects [2^, 21 1: 
For instance, Ref. [21j supposed the existence of large 
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scale structures in the GMF aligned with the sources, 
leading to magnetic funneling of CRs. Anisotropic MHD 
turbulence in the interstellar magnetic field has also 
been proposed as a candidate to create such small scale 
anisotropics in the CR arrival directions at Earth (2^ . 
On the contrary, Ref. [2^ assumed that CRs may be ac- 
celerated very locally from magnetic reconnection in the 
magnetotail. 

Small scale anisotropies from local magnetic turbu- 
lence. — In this Letter we show that energy-dependent 
medium and small scale anisotropies necessarily appear 
on the sky, provided a large scale anisotropy exists, ei- 
ther dipole -for instance from the inhomogeneous source 
distribution- or dipole and quadrupole, such as in Ice- 
Cube data [3] . The small scale anisotropies are due to the 
structure of the local turbulent GMF, typically within the 
scattering length X{p) from Earth. Since the turbulent 
field modes relevant for cosmic ray scattering depend on 
the CR rigidity and, therefore, the local volume of the rel- 
evant turbulent field is energy- dependent, the small scale 
anisotropies must be energy- dependent. Any -isotropic- 
turbulent magnetic field would generically produce small 
scale anisotropies. The diffusion approximation cannot 
describe them because it averages over different turbulent 
magnetic field realizations. It thus averages out small 
scale anisotropies that are uniquely created by the con- 
crete local realization of the field. Neither anisotropic 
turbulence, nor directional correlations with source di- 
rections on the sky, are required to explain the data. 

CR trajectories in three-dimensional turbulent fields 
can be regarded as random walks. For times » X[p)/cq 
the CR transport can be considered as a Markovian pro- 
cess where CRs diffuse on random scattering centers. Let 
us consider a point-like source emitting particles in all 
directions. For distances from the source < A(j3), the 
propagation in the local turbulent magnetic field still has 
memory: The particle trajectories are locally determined 
by their initial directions, and a very small change of the 
initial angle would not lead to very significantly different 
trajectories for distances < \{p). For distances to the 
source > A(p), the process becomes Markovian: CRs on 
nearby trajectories virtually lose memory of times longer 
ago than ~ X{p)/co, and their propagation is uncorre- 
lated with initial conditions. 

Since Galactic CR energy loss can be neglected and 
CRs can be back-tracked, the situation is the same for a 
point-like observer of size ^ such as Earth. There- 
fore, the small scale anisotropies that are observed at 
Earth arc due to the propagation within the very local in- 
terstellar turbulent magnetic field, within the "sphere" of 
radius ^ X{p). The hot spots at Earth arc regions where 
particles are statistically more connected to parts of the 
"sphere" where the dipolar flux impinging from outside 
the sphere is larger, and cold spots are connected with 
the part of the external dipole which has a deficit. We 
also note that the motion of the Sun in the Galaxy cannot 
smooth out the anisotropies because the distance traveled 



during the experiment lifetime is ^ A(TeV). Moreover, 
the local turbulence configuration cannot change notice- 
ably on time scales of ^ X{p)/co. 

To put this more quantitatively, let us imagine two 
spheres around Earth with radii r ~ A(p) and i? 3> A(p), 
respectively. Back-tracking trajectories of total momen- 
tum p from Earth in a given direction parametrized by 
the unit vector n will give a unique point on the sphere 
of radius r, represented by the unit vector Gp(n). The 
function Gp is determined by the magnetic field within 
the small sphere and depends on p. It is smooth because 
trajectories are still non-Mar kovian on scales < X{p). 
When back-tracking further, memory of initial conditions 
is largely lost, in particular when n is averaged over the 
experimental angular resolution. Still, the expectation 
value of the crossing point i?nfl after diffusion within 
the large sphere must be (Rnn) ~ A(p)Gp(n). There- 
fore, the probability pp(n, rii?) for a particle emitted in 
direction n at Earth and back-tracked to cross the large 
sphere at Rnji can be estimated as 



1 + Gp(n) • ur 
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which is normalized, j pp{n,nii)(Pnii ~ 1, when inte- 
grated over the sphere. If Fii{n]i) is the differential in- 
tensity impinging from outside the large sphere at point 
Rnji, the flux F{n) seen at Earth in direction n will then 
be given by 



F{n) = j FR{nR)pp{n,xiR)(fnR. 



(6) 



We here assume a uniform gradient of CR density such 
that Fr is the sum of an isotropic piece _Fo ^iid a dipole 
of relative strength A in z— direction e^. 



FR{nR) = Fo [1 + A [nR ■ e,)] 



(7) 



This corresponds to 'Wn/n = (A/R) and, according to 
Eq. ([2]), results in a dipole d{p) ~ A(p)(A/i?) observed 
at Earth. At the same time, the flux at Earth, Eq. ([6]), 
contains multipoles aim = J F{n)Yim{n)(fn/Fo, nor- 
malized to Fq, given by 



(8) 



aimip) = |^(p)| / Gp(n) • Y'/„i(n)d n. 



where Yimin) are the general spherical harmonic func- 
tions. Eq. ([S]) follows after performing the integration 
over ur and expressing A in terms of the dipole ob- 
served at Earth. In the absence of deflection within 
a distance < X{p) from Earth, Gp(n) = n and with 
yio(n) = [3/(4^)]'/' n • e„ Eq. ® gives ai^ip) - 
<^/i^™o|^(p)| [(47r)/3]'/'. This yields F{n) = |^(p)|i^on- 
e^, corresponding to diffusion with only the dipole 
present. As in a gas, particles may need several scat- 
terings before losing memory of their momentum. Since 
the positions of scattering centers in gases are time- 
dependent, the time-averaged behavior as seen by a given 
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FIG. 1. Renormalized CR flux predicted at Earth for a 
concrete realization of the turbulent magnetic field, after 
subtracting the dipole and smoothing on 20° radius circles. 
Primaries with rigidities p/Z = 10^® eV {left panel) and 
5 X 10^^ eV {right panel). See text for the field parameters 
and boundary conditions on the sphere of radius R — 250 pc. 



observer is equivalent to particles traveling on straight 
lines of length ~ X{p) {i.e. Gp(n) = n) and losing 
memory after each scattering, also known as "molecu- 
lar chaos" . In the case of a concrete local turbulent field 
realization within ^ A(j>), the function Gp(n) ^ n in 
general has structures on all angular scales. Eq. (|8]) pre- 
dicts \aim{p)\ ^ \^{p)\' Clearly, higher multipoles should 
be of the same order as the dipole impinging from outside 
the sphere around Earth. As shown in Fig.[TJ the smaller 
scale anisotropies are much more rigidity-dependent than 
the dipole. Since the CR composition is not pure proton, 
and since experiments present data with "broad" energy 
distributions and uncertainties Ap/p ~ 0.3, small scale 
anisotropies sum up in a non-constructive way, which 
typically leads to smaller amplitudes than that of the 
dipole by a factor of a few to ten. If the structure of 
the magnetic field and thus the function Gp were known 
within a distance X{p) around Earth, Eq. (jH]) would 
even allow to make concrete predictions for the small 
scale anisotropies. It is also clear that the predicted 
anisotropies should not depend significantly on i? > X{p) 
and should thus converge quickly with increasing R. Be- 
low we confirm this with numerical simulations. In prac- 
tice, the Alfvenic part of the turbulence is anisotropic 12411 
and scattering is mostly dominated by fast modes [25| . 
The condition Gp(n) ^ n always holds, which systemat- 
ically ensures the appearance of small scale anisotropies 
even if scattering is not isotropic. If the diffusion tensor 
averages to isotropic diffusion on scales ^ X{p), Eqs. ^ 
ff. are unchanged. 

Numerical simulations. — We back-track CRs in con- 
crete turbulent magnetic field configurations. The tur- 
bulent magnetic field is generated on three-dimensional 
grids with the method presented in Refs. [H-il. We 
use a Kolmogorov spectrum between L,nin rg and 
imax = 150 pc and with an rms strength -Brms — 
CRs are injected with isotropically distributed random 
directions from Earth, located in the Galactic plane 
XY at (0, 8.5kpc), until they cross a sphere of radius 
R = 250 pc > imax around Earth. Trajectories are con- 
sidered as fiux tubes and a CR crossing the sphere at po- 
sition RnR is weighted by Eq. d?]), except that here the 
dipole is taken towards ey . We take | Vn/n[ = 1/(290 pc) 



&ndp/Z = 10, 50PeV, so as to detect the predicted effect 
with reasonable statistics. 

For this set of parameters we find the expected dipole 
strength \6{p)\ ^ 6% and direction at Earth. Its exact 
direction and amplitude slightly depend on the concrete 
realization of the local turbulent field and on the position 
of the observer, for the same reason as the appearance 
of small scale anisotropies. Since the dipole scales as 
\5{p)\ oc |Vn/n|p^/'^ for Kolmogorov turbulence, taking 
a more realistic gradient |Vn/n| ~ 1/(1 kpc) and energy 
p ~ lOTeV would rescale its amplitude to ~ 0.1%. We 
subtract the dipole from the predicted fiux at Earth and 
show the relative residual intensity maps in Fig. [1] after 
smoothing on 20° radius circles on the sky. The statisti- 
cal fluctuations in these computations are below ~ 0.5%. 
One can see statistically significant features of various 
amplitudes and shapes, some of which may well resem- 
ble the data. Their significance trivially depends on the 
smoothing radius. Varying the magnetic field realization 
we find that, while the dipole only slightly varies, the 
shapes and positions on the sky of the small scale features 
are strongly realization dependent. In the left panel of 
Fig.[T]there is a "hot spot" with twisted shape and ampli- 
tude ~ 6% in the lower left quadrant and a "cold spot" 
(~ —6%) in the upper right quadrant. As predicted, 
the amplitude is comparable to that of the dipole. At 
five times larger rigidity the features are strongly differ- 
ent, as seen in the right panel of Fig. [TJ For instance, 
the cold spot at lOPeV in the upper right quadrant is 
transformed into a hot spot. The energy-dependence of 
features may also account for the different spectrum [1] 
seen in the Milagro hot spot. The amplitudes of fluctu- 
ations at 50 PeV are larger than at 10 PeV here, because 
we kept the same gradient. 

Fig. [2] shows how the predicted anisotropies in Fig. [1] 
depend on the maximal back-tracking distance R. At 
10 PeV and 50 PeV, convergence of the sky maps is es- 
sentially achieved for i? > 25 pc and R > 50 pc, respec- 
tively. These length scales roughly correspond to the val- 
ues of X{p) in the vicinity of Earth in the given magnetic 
field realization. Averaged over many realizations, for 
6 ~ 1/3, their ratio, as well as the ratio of corresponding 
anisotropy amplitudes in Fig.[Tl should equal 5^/^ accord- 
ing to Eq. ([3|), and Eqs. (g]) and dH), respectively. The 
small deviation from this average of the ratios simulated 
in the given realization is thus due to "cosmic variance" . 
Fig. [2] shows that the small scale fluctuations arise from 
the local field within ^ X{p), as predicted above. 

Fig. |3] (first panel) presents CR flux anisotropies in a 
30° X 30° sky patch, after smoothing on 5° radii circles. 
The three other panels show the trajectories of four CRs 
arriving at the red crosses in the flrst panel (two cho- 
sen in an excess region and two in a deflcit region). The 
third and fourth panels show that, at distances i? ^ A(p) 
from Earth, the two trajectories arriving in the hot spot 
tend to come from the direction of the CR density gra- 
dient {Y > 8.5 kpc) while the other two come from the 
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FIG. 2. Same as Fig. [T] for boundary conditions imposed on concentric spheres around Earth with radii R — 100, 50, 25, 10 pc 
(resp. first, second, third and fourth columns). Upper row: p/Z = 10^^ eV; Lower row: p/Z = 5x 10^'^ eV. 
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FIG. 3. First panel: Renormalized CR flux predicted at Earth for p/Z = 10^^ eV in the sky patch / = 150° - 180° and 
6 = 0° ~ 30°, smoothed on 5° radius circles. This is a blow-up of Fig. [1] (left panel); Second, third and fourth panels: 
Trajectories of four CRs back-tracked from Earth, projected onto the XY plane and within ~ 200, ~ 20, and a few pc 
from Earth {resp. second, third and fourth panels). The CR initial directions, denoted by red crosses on the first panel, are 
I — 154°, 156°, 163°, 164° (see values in panel labels) and 6 = 5° in celestial coordinates. 



opposite direction, consistent with Eq. ([5]). On larger 
scales (second panel), the initial directions arc more uni- 
formly distributed, again consistent with Eq. ([5]): The 
small scales reflect the last part of the particle trajecto- 
ries (< A(p)/co), before they are detected by the point- 
like observer. 

We verified that small scale anisotropics also appear 
for anisotropic CR scattering, by performing simulations 
with an additional large scale field. 

Conclusions. — We have shown that the observed inter- 
mediate and small scale anisotropics in the Galactic CR 
arrival directions can be naturally explained as the conse- 
quence of CR propagation in a turbulent magnetic field. 
The observed anisotropics could thus be one of the first 
direct manifestations of the turbulent Galactic magnetic 
field within the scattering lengths of TeV-PeV CRs, and 
thus within a few tens of parsecs from Earth. Formally, 
this effect could have similarities with the "CR scintilla- 
tions" in the inner heliosphere [1^. In the future, this 
should allow new insights into the CR transport in our 
Galaxy and contribute to our knowledge of the structure 
of local -and notably interstellar- magnetic fields. 
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